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In olie preGerni vro3.'k, 
dis-uribivuioiis in Pluii,i,e~ciru 
B ini '’go iD-omcni iieij^od, Bor 


temperatures and residual stress 
grinding lias oeeii obtained using 
simplicity, a tifo dimensional 


approximation for the grinding process is made. Transient 
heat conduction equation was solved to obtain the temperature 
distribution in the workpiece. The motion of the grinding 
vrheol has been simulated by a step-by-step movement of the 
heat jiluGc over the wheol-worl:pi.ece contact cone. The above 
t emiierature historjr ''•i.as been used to find the residual stress- 
es in the subsequent elastic-plastic aialj^sis , 

The numericoh results obtained, confirms that the 


tensile residual streosos may reach the ultimate strength of 
the workpiece material. This causes grinding cracks whose 
presence reduces the fatigue strength of the ground part 
resulting in its pr-miEiture failure. An increase in the depth 
of cut and/or the table speed results in increase of the peak 
value of the tensile residual stresses, Variation in heat 
transfer coefficient, so as to simulate the effect of cooling 
fluid, has marginal effect on the peak valut3 of the tensil; 
residual stress. The results so obtained compare well with 
the experimental data available in the literature. 
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IRTRODUCTIOIT MD LIPSRAIDUilE REVIEW 

1.1 INIRODUG'IIOR 

G-rinding is a metal cutting process which is similar 
in many ways to other commonly employed methods of metal 
removal such as turning, milling or shaping* In grinding, 
material removal takes place due to shearing just as in other 
cutting operations. Ihe grinding wheel may be described as a 
multi-toothed milling cutter, each tooth consisting of a small 
abrasive particle. Of all the metal cutting processes, grind- 
ing is undoubtedly the least understood. Because of the random 
grit geometry, high cutting speed, and small depth of cut 
taken in grinding, mechanisms of grinding are diffici0.t to 
observe and evaluate. 

The premature failure, in service of components manu- 
factured by the grinding process and the appearance of surface 
cracks during grinding of certain materials have made it an 
economic necessity to carry out detailed investigations into 
the surface integrity aspect of the process, including resi- 
dual stresses. When a grain first contacts the worlspiece 
surface (Pig. 1.1), sliding and ploughing occurs. As the 
grain moves it experiences more and more uncut material and 
the forces required to remove this material increase, and 
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ultimately the material is removed in the form of a chip. 
During this process of material removal, large amount of 
heat is generated in the workpiece, due to plastic deforma- 
tion of the work material, due to frictional sliding of the 
grain over the workpiece surface, and due to the rubbing 
action between the chip and the grain. Most of the heat 
generated (about 80-85% ) enters the worl^piece and the 
resulting temperature rise may cause burn or thermal cracks 
on the workpiece surface. The above process of material 
removal leads to inhomogeneous plastic deformation. Therefore, 
individual elements of the workpiece undergo different 
changes in their dimensions. Besides the forces in the work- 
piece, which are in equilibrium with the applied forces, 
additional forces 'maintain equilibrium within the workpiece 
after the release of external loadings. These internal forces 
cause stresses which are called residual stresses. 

The residual stresses can have three different 
origins : 

(l) Borces between the worlspiece and the tool 

The forces in a cutting process cause elastic and 
plastic deformation. The elastic deformation subsides as 
soon as the cutting process terminates. The plastic defor- 
mation remains and results in compressive residual stresses 
in the workpiece surface. 



(2) Heat Expansion 


Ihe material expands locally owing to the heat 
generated dijring the cutting process. If the cutting tempera~ 
tures are too loigh, inhomogeneous expansion will exceed the 
limits of elastic deformation of the material and leads to 
tensile residual stresses after cooling. In extreme cases, 
cracks can appear on the workpiece surface* 

(3) Structural Transformation 

G-rinding temperature, sometimes reaches the hardening 
temperature of the material. Because of the temperature 
gradient, tempering occurs in certain layers below the surface, 
while the workpiece surface is rehardened by virtue of reaching 
the critical cooling speed throu^ the application of coolant 
liquid. The tempered material has a smaller volume compared 
with the original lattice, whereas the newly-hardened workpiece 
material reverts to same volume as before. Tensile residual 
stresses are thereby created. 

Under fine grinding conditions, the temperatures 
developed are well below the temperature that leads to struc- 
tural transformation. The possibility of residual stresses 
due to this transformation is therefore negligible. The 
effect of mechanical loading on residual stresses were studied 
by Mishra et al, [17j, and it is found to be negligible. 
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Residual stresses are therefore mainly due to non-uniform 
temperature distribution in the worlspiece, generated during 
grinding. 

fhe thermal gradients that are imposed on the upper- 
most 3.ayers of the body would result in a large thermal 
expansion. Becuase of the deformation constraints created 
by the much stiffen and not thermally affected material of 
the workpiece, large compressive stresses are built up during 
heatlng-up period, which may surpass the yield limit of the 
material, especially at elevated temperatures. During cool- 
ing down period, the upper layer tends to contract elastically 
more than the substrate, creating tensile residual stresses 
in the uppermost layer. 

Residual stresses are detrimental to the fatigue life 
of the ground part, previous studies have indicated that 
compressive residual stresses increase the fatigue strength 
of ground parts. Since grinding operation is the last process 
to be performed over the part, the stresses that are created 
during this process determines the life of the part, 

1.2 lITBRilURE SURVEY 

URIEDING SBMPERilURES 


The first attoupt to estimate the grinding temperatures 
appears to have been made by Guest [1]. He proceeded with the 
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assumption that half of the total heat generated passes away 
with the chips and related the specific horsepower to the 
temperatures. But due to the complexity of the process and 
difficulties involved in measuring the temperatures accura- 
tely, analytical approadies have Been resorted to, These 
temperatures models [2-4] employed the moving heat source 
theory proposed hy Jaeger [5]« Shaw et al, [2] used this 
theory, with the assumption, that an average grit on the 
wheel surface has a raice angle of aero degree, and estimated 
the worlipiecG surface temperatures for the fine surface 
grinding. They also measured the temperatures experimentally 
and concluded that as table speed increases, temperature 
increases, Sato [6], determined the temperature at the wheel 
work contact surface from the heat quantity conducted into 
the worlspiece, which he found out analytically. His theory 
was based on the assumption that the workpiece is a semi- 
infinite body and the grinding wheel is an instantaneous heat 
source. He also found out that 85% of the grinding energy 
is conducted as heat into the workpiece, Dos Ruisseaux and 
Zerkle [5] calculated the grinding temperatures as the sum 
of a local temperature due to grinding by an individual grain 
and a grinding zone temperature due to grinding by all other 
grains. The relations were established by considering the 
grinding geometry. Malkin and Anderson [7] considered the 
total grinding energy as composed of energy required for chip 
formation, ploughing and sliding. He assimed that a portion 
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of this energy is conducted as heat into the -workpiece and 
related these results to grinding temperatures and workpiece 
hurn, Sauer [8], in his theory combined both experimental 
and analytical techniques, and calculated the temperatures. 

Ihe effect of coolant and rise in temperature due to plastic 
deformation was considered. 

RESlDUil STRESSES 

Several researchers have attempted to determine 
residual stresses experimentally. A considerable portion 
of the work is attributed to Gilkman and Stepnov [9], Frisch 
and Thomsen [10], Marshell and Shaw [ll] , Letner [12,13], 
Colwel, Sinnot and Tobion [14], Halverstadt [15] and others. 

A review of existing experimental work brings out two main 
aspects. Firstly the experimental determination is often 
done by bending-deflection method. This method involves 
removing of material from the stressed surface electro- 
chemically or by etching, whilst measuring the curvatures at • 
regular intervals using a sensitive measuring device. This 
method is beset with difficulties in controlling the unifori 
mity of layer removed, and accurate measurement of curvature 
of the specimen, which effects the stresses calculated. 
Secondly, the stress distribution is influenced by a number of 
grinding parameters such as wheel speed, wheel grade, depth 
of cut etc., besides the met all in* gl cal nature of the material 
being ground. Both tensile and compressive residual stresses 
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have hccn found near the surface under differ oat sets of 
conditions, and published data are often contradictory, 

Thus it is difficult to fona a clear picture as to the 
influence of processing parameters on the residual stresses 
introduced. 

Other methods that have been employed for estimating 
the residual stresses are X-ray diffraction method, Sachs 
boring method etc. [16], The measurement of residual stress- 
es by X-ray, utilises the interatomic spacing of surface 
lattice planes as the guage length for measuring strain. 

Using Bragg’s law, which expresses the relationship between 
the distance from each other of a given set of lattice planes, 
the wavelength of X-ray diffraction, the order of diffraction, 
and the measured diffraction angle, the usual stress equation 
for the two-exposure technique is given by 


where 
A© = 


K 


3 


E Aq. 

^ ^ sin^ ¥ tan 6 


KA © 


change in the Bragg an^e between normal and inclined 

exposures 

constant 
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This method is particularly useful, where surface stresses 
are important, and since the exposure area is small, steep 
stress gradients can he obtained. But the disadvantage is 
that the equipment is expensive, and if the area to be exa- 
mined is large, then measurements are tedious. This method 
is inapplicable to specimens that have been subjected to 
uniaxial plastic deformation and difficult to apply to heat- 
treated high-tensile steels and severely cold worked materials. 

Sach’ s boring method employs the technique of boring 
out a cylinder or tube in stages and measuring longitudinal 
and circumferential strains at the outer surface caused by 
the release of residual stresses. The assumptions involved 
in this method are, material is isotropic and Youngs modulus, 
Poisson’s ratio ai-e constant, residual stresses are distri- 
buted with rotational symmetry about the axes, the tube formed 
by boring is circular in section and the inner and outerwall 
surfaces are concentric, specimen is sufficiently long to 
prevent lateral bending. The disadvantage over the bending 
deflection technique is that plastic deformation should not 
occur on layer removal, strain guages used must be hi^ly 
sensitive. 

The residual stress distribution has been obtained 
analytically by Mshra et al. [17] using Finite Element Method. 
They have obtained the temperature distribution in the body, 
from the r Ration given by Jaeger [5]» and proceeded to obtain 
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stresses using elasto'-plastic analysi.s. No attempt -was 
made to evaluate the effect of grinding variables on resi- 
dual stresses, and comparison with the experimental values, 

lo3 FINITE BLli-ISIT TECHNIQUE 

In many practical situations, the geometry of the 
S 3 ’'stem and the boundary conditions are too complex to yield 
analytical solutions <, Such problems can, however be solved 
by numerical methods. Popular amongst these numerical 
techniques are Finite Difference and Finite Element Methods, 
The Finite Difference scheme of a problem gives a pointwise 
approximation to the governing equations. The accuracy of 
this model can be increased by taking more number of grid 
points, but its main drawback is, it becomes hard to apply 
this technique for problems with irregular geometries or 
unusual specifications of the boimdary conditions. Whereas 
Finite Element Method envisions the solution region as built 
up of many small interconnected sub-regions or elements. 

This method gives a piecewise approximation to the governing 
equations. The basic premise of the Finite Element Method 
is that a solution region can be analytically modelled or 
approximated by ri^lacing it with an assemblage of discrete 
elements. Since these elements can be put together in a 
variety of ways, they can be used to represent exceedingly 
complex shapes. 
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111 more and more engineering situations today, it 
is necessary to obtain approximate solutions to the problems 
rather than exact closed form solutions, Ilic Finite Element 
Method is a numerical analysis technique for obtaining 
approximate solutions to a wide variety of Engineering problems, 

I.! a conti’na;-.'! problem of any dimension, the field 
variable (whether it is pressure, temperature, displacement, 
stress, or some other quantity) possesses infinitely many 
values because it is a function of each generic point in 
the body or solution region. Consequently, this problem is 
one with an infinite number of unknowns ; The finite element 
discretization procedure reduce the problem to one of a 
finite number of unknowns by dividing the solution region 
into elements and by expressing the unknown field variables 
in terms of assumed approxiE.ating functions within each 
element. The approximating functions (sometimes called 
interpolation functions) are defined in terms of the values 
of the field variables at specified points called nodes or 
nodal points. The nodal values of field variable and the 
interpolation functions for the elements completely define 
the behaviour of the field variable witliin the elements. 

The nature of the solution and the degree of approximation 
depends not only on the size and the number of elements 
used, but also on the interpolation function selected, A 
Finite Element Model of a problem gives a piecewise approxi- 
mation to the governing equations. The important feature 
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of the Pinite Elemont Method is its ability to formulate 
solutions for individual elements before putting them together 
to r (present the entire problem, 

The various steps involved in the solution of a 
problem using Finite Element Method are : 

(1) Discretization of continuum, 

(2) Selecting interpolation functions, 

(3) Finding the element properties, 

(4) Assembling the element properties to obtain the 
sys t em equ ati ons , 

(5) Solve the system equations, 

(6) MaFe additional computations for stress/strain 
evaluation, 

1,4 PRESENT wore: 

A frequent cause of failure of ground surfaces is a 
network of practically invisible cracks that develop during 
or immediately after grinding as a result of considerable 
heat generated during the process, and the presence of hi^ 
residual stresses in the ground surface. If the tensile 
residual stresses produced by grinding are high enou^ to 
exceed the rupture strength of the material, cracks -will form 
under the operational loads to which the surface is subjected, 
this crack propagates leading ultimately to the failure of 
the part. 



In the present work, an attempt has been made to 
study the effect of grinding parameters on the temperatures 
generated and the residual stresses thus formed. Por 
simT>licity, a two dimensional case for plunge cut grinding 
has been considered, temperatures were obtained by solving 
the basic heat conduction equation in a solid using Finite 
Element Method, The cooling effect is introduced as a speci- 
fic thermal boundary condition at the workpiece surface. 

An iterative procedure is employed for the step-by-step 
movement of the temperature field, in order to simulate the 
movement of the grinding wheel over the workpiece surface. 
Since the residual stresses are mainly duo to non-uniform 
temperature distribution in the workpiece an attempt has been 
made to obtain the clastic-plastic stress field within the 
workpiece from the temperatures calculated above. The results 
obtained using this approach have been compared with the 
available experimental data. 






1,4 


.Jt'i ' ^ 


- 4 . 






Fig, 1,3 Two dimensional thermol mode! representing 
grinding action . 




CHiPTBR II 


K-IERll&L Al'JiiYSiS 


2,1 IFIRODUCIIOn 

A large percentage of energy used in grinding is 
transformed into thermal energy ■which causes the workpiece 
temperature to increase* For obtaining better grinding 
resTolts, it is necessary to minimise this temperature rise 
because of the foliating % 

(1) The maximum temperature in the contact area 
between the wheel and the workpiece determines the rate of 
jhemical reaction at the interface. 

(2) The maximum temperature of a finished wor!kpiece 
surface determines the integrity (or micro structural damage) 

)f the surface. 

(5) The a'verage temperature of the bulk of the work— 
liece determines the dimensional accuracy after grinding* 

Most of the theoretical models for obtaining grinding 
lemperatures are either based on Jaeger^ s band source model 
3 r moving line source model [5]. Jaeger’s thermal model, 

1 band heat source is assumed to move on the workpiece surface 
it a speed equal to that of the worlspiece velocity v. Bas,^ on 
ihis model Des Ruisseaux and Zerkle [3] calculated the surface 
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■bemperatures ’by superimposing tho temperatures due to a 
single grain and due to all other grains in the contact zone, 
Sauer's analysis [8] also considered the heat input due to 
plastic deformation of the workpiece hut this heat due to 
pla^stic deformation in often of secondary importance. Sato’s 
theory [6] was based on the temperature rise of the surface 
of a semi-infinite body under an instantaneous heat source. 

He calculated the heat generated during grinding from the 
tangential grinding forces. 

2,2 iNllYSIS 

The temperature model is characterised by three sets . 
of fundamental physical quantities, (l) The heat input into 
the worlspiece and its distribution over the contact area, 

(2) the real contact length (1 ) and (3) the thermal charac- 

G 

tenistics of metal (conductivity k, specific heat per unit 
volume ^c) • 

(l) Actual heat input 

The fraction of the total energy dissipated by the 
heat source flowing directly into the workpiece is of special 
importance from the point of view of thermal load on the work- 
piece, This was experimentally measured by several authors 
[6,7,8], Roughly 70-90% of the total consumed energy flows 
instantaneously into the workpiece causing a rapid increase 
in the local workpiece temperature. A large specific metal 
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removal rate usually reduces this percentage hut blunt wheels 
show a tendency of increasing tliis value [7]. In the present 
work, it is assumed that 80X of the total energy dlfcssipated by 
the heat source flows into the workpiece instantaneously, 

(2) Contact length 

The real contact length 1 is a function of wheel 
dimensions, composition and characteristics. It also depends 
upon the kinematic parameters (speeds and forces) because they 
govern the normal grinding force component and elastic flatten- 
ing of the wheel in the contact area. The length 1 is approxi- 
mat el y equal to 



( 2 . 1 ) 


where 

D = effective wheel diameter (mm) 
d = depth of cut (mm) 

The width of the interference zone is equal to the width of 
contact between the vrheel and the workpiece. For plunge 
cut grinding it is equal to the width of the workpiece. (Pig. 
2.1),* Heat flux distribution in the grinding wheel-work 
interface has been studied by Jaeger [5] and Sauer [8]. They 
found that this flux distribution has effect upon the maximum 
temperatures developed during grinding. Further, only a few 
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microns below the groimd surface arc most affected by the 
shape of the flux distribution. So a uniform heat distribution 
may be assumed for the evaluation of temperatures without 
affecting very much the overall accuracy of the results. 

(3) Thermal characteristics 

The thermal characteristics of workpiece material, 
in general, are temperature dependent. However, an adequately 
chosen average value may give a fairly good approximation of 
the temperature distribution as reported by Sauer [8], This 
assumption simplifies the solution of the heat transfer equa- 
tion for evaluating the temperatures using Finite Element 
Model. 

The temperature rise in the workpiece can be calculated 
by approximating the combined action of all grains and by 
considering the wheel as a single heat source moving with a 
velocity, V, over the workpiece. The bulk temperature rise in 
the workpiece due to one pass of grinding wheel, assumed 
initially to be at uniform temperature, T^, is obtained from 
the temperature rise due to interference zone surface heat 
source* 

The interference zone boat flux is given by [8] 


Q- 


•F 


a F Y 

rx> 


t 


( 2 . 2 ) 
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vrliere 

a = fraction of the total grinding energy entering into 
the •workpiece 

= horizontal grinding force (grnf) 

Y = wheel speed (mm/sec) 

B = ■frri.dth of the T/rorlq^iece (ram) 

J = JoifLes constaait ( gmf-mm/cal ) 

2- o 

Hewtonian cooling at a uniform rate h (cal/mm -sec- C) from 
the surface is assumed for calculating hulk temperatures. 

Two dimensional transient heat conduction eauation, 
for a uniformly tiiick worlq)iece of thickness t is given hy [18] 


1 _ 

3x 







( 2 . 3 ) 


The desired temperature must satisfy the following Boundary 
conditions i 

(a) T = TCxjyjv) on S^, i: > o (Pig. 2.2) 

Co) ISc M 'ly + ^ + h(T-T^) = 0 on Sj, t > o 

(c) T = T^Cxjy) in D, i; = o 


where 

k ,k = theraal conductivities in x and y directions, respecti- 

jS. ^ 

veiy ( cal/iBii>-S0C— C) 


(V 
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Q = internal heat generation 

I'q = initial temperature in the domain (°C) 

= specified boimdary temperature distribution (°C) 
= direction cosines of the outward drawn normal 
vector n to the bounding curve 
Op = heat flux lost or added at the boundary due to 

p 

conduct! on ( cal/ sec~mm ) 

11 ( 1 - 1 '^ ) = heat lost at the boundary due to convection to 

ambient temperature, T , with convective heat 

00 

transfer coefficient h 

p = density of the worl-cpiece material (gm/mm^) 
c = heat capacity of the material (cal/gm~°C) 

Equation (2,3) in its variational form can bo written as 


t 



+ y Op 1 dEdZ + l)dBdZ . 

o 0 B 

( 2 , 4 ) 

The entire continum is discretized into a finite number of 
sub-regions as shown in Fig, 2,2, Minimisation of this func- 
tional requires an explicit formulation of temperature T in 
each element. 

Assuming the temperature distribution within each 
element as 


^ 4 
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r 

I N^(x,y)T^(T) 
i=l 


[U(x,y)] {I(t)}(®^ 
(2.5) 


where 

r = number of nodes assigned to each element, 

~ discrete nodal temperatures, 

iU2^jy)= shape functions or interpolation function., 

The functional has a minimum, when its first variation with 
respect to T vanishes* The resulting equations in the matrix 
form can be written as [19] 

/ / [IS: {“.xJ H,y}l>I,y}]tady + 

gCe) 

f f ”0 {N} [H] dxdy{ = {Q}^®h 

s(-) 

- [^1 (I) + {k ly (®^ , (2.6) 

where superscript e refers to the element considered and S' 

1 „ are the partial derivatives of the shape function or inter 

»y 

polation function. After subsequent integration over each 
element and assembling all the element matrices, to form the 
global matrix, leads to a system of first order differ^tial 
equation. This can be represented in matrix form as 
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[Kq] {i} + {I } = { Q) - _ [Kg]{ S}+{EtJ, 

(2.7) 

where 

= thermal capacitance matrix, 

[K^p] = thermal conductivity matrix, 

{Q} = internal heat generation, 

^ ~ heat loss or gain due to conduction, 

[%] = f h{Hg}[iyaB 

1 Heat loss due to convection 

[K^ ] = f h T„{H } dB boundary, 

S 2 

Hg = interpolation functions, on the houndary. 

The boundary conditions that were considered in the present 
case are, cooling at a uniform rate by convection with convec- 
tive coefficient h ( cal/mm^-sec-°C) from the surfaces AB, BC 
and AD (Big, 2,5) • ' The surface DC is insulated so that no heat 
can pass through the workpiece into the magnetic chuck on which 
the workpiece is mounted. 

The boundary conditions are th®a substituted and the 
resulting first order differential equations are solved using 
Grank-ilich Olson recurrence relation giving [20] 

( 2 . 8 ) 
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xfhcre P is the R.II.S, of the equation (2,7)« Here it is 
assumed that the temperature varies linearly with time 
durin£; the arhi^^ary time interval At, The time intervalA 
is obtained in the following manner. 

The movement of grinding wheel is approximated such 
that, it moves in discrete steps* The time for which the 
interference zone heat flux is acting on the workpiece, i.e., 
the time required for the wheel to travel between two success- 
ive steps is 



For this time period, it is assumed that the heat flux is 

stationary. This time period is suitably divided into a 

number of time intervals At, and equation (2.8) is solved. 

After the time period t ^^ ^ , the wheel is shifted to a. distance 

1 (i.G. the heat flux zone is shifted) and then equation (2,8) 

c 

is again solved for a period t ^. using the temperature 
obtained in the preceding case. This process is continued 
till the contact between the wheel and the workpiece ceases. 
Thus the temperature history in the workpiece is obtained, as 
the wheal proceeds. 

2.3 COMPARISON ¥ITII EXPERIMENTAL RESULTS 

A computer programme was developed to solve equation 
(2*7) by the recurrence relation given in equation (2.8), and 
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the grinding temperatures were obtained for different sets 
of conditionSo fhe model was tested by comparing with the 
experimental results obtained by Sauer [8], A close agree- 
mont between the theoretical and experimental results is 
obtained. A typical temperature distribution curve is shown 
in Pig. 2.15. Ihe input data used is given in Table 2.1A. 

2.4 numiCAI RESUITS 

For the present work, the data was taken from the 
experimental work of Pande and lal [21] and is tabulated 
in Table 2. IB. The resulting temperature distributions, 
are shown in Pigs. 2.4 to 2.14. Prom Pig. 2.4, it can be 
seen that the maximum temperature attains a constant value, 
and as the wheel moves over, a wake of temperature field is 
left belli nd. The temperature distribution in the contact 
length (1.) shows a steep increase to the peak value and 
then drops rapidly, confirming the theoretical predictions 
of Jaeger [5]. The peak temperature shifts to the rear of 
the heat source Icngtli, as the contact length increases. 
Tiiis is also in accordance with Jaeger’s model. 

The temperature distribution inside the workpiece 
at various depths from the surface is shown in Pig. 2.6. 

As the d^th is increasing, the peak temperature shifts 
towards the rear side of the heat source length. Pig. 2.7 
to 2.14 show the plots of temperat'ure distribution for 
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vai'ious dci^tli of cuts and talDle speeds. As the table 
speed inci’casos, there appears to bo a slight increase 
in the peak teroperaturos, lig. 2.17 > tut at lower 
table speeds tho depth of heat affected zone increases. 
I'liis is shown in Figs, 2,6 to 2.8, 

She effect of depth of cut is more pronounced on 
the temperatures. As the depth of cut increases, the peak 
temperatures increase significantly. These are shown in 
Figs, 2.8-2,10. 
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Fig. 2.7 Temperature distribution at various depths 
from the surface (set no. 2). 
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( a ) y » 0 mm 

(b) y » 0.3 mm 

(c) y * 1.0 mm 

(d) y » 1.5 mm 



Dirtction of whool movgment , mm 


Fig. 2.8 Temperature distribution at various 
depths from the surface, (set no. 1). 
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900 (q) y ^Offtfin 



Direction of wheel movement , mm 


Fig. 2.9 Temperature distribution at various depths 
from the surface (set no. 4). 
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Fig. 2.11 T«mp«rature distribution nt various depths from 
the surface (set no. 8). 
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Fig. 2.13 Ttmpiraturt distribution at various depths from 
the surface (set no. 5). 
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(a ) y * 0 mm 



Fig. 2.14 Temperature distribution ot various depths- 
from the surface (set no, 6). 
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Effect of depth of cut on surface 
temperature. 



Table speed, m/min 

F:g. 2.17 Effect of table speed on surface 
temperature. 
















TABLE 2.1A I lOTUT DATA EOR CALCULATING TEi€> MATURES. COI^ ARISON WITH SAUER 
EXPERIMENTAL RESULTS [8] 
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CEAPTBR III 


RESIDUAL STRESS MADYSIS 

3.1 INTRODUCTION 

As explained earlier, the residtial stresses are 
mainly due to non-uniform temperature distribution in the 
wortpiece developed during grinding. These stresses, 
which are tensile in nature, may sometimes surpass the 
yield limit of the workpiece material, resulting in crack 
formation. These cracks, however small they may be 
initially, grow under the operating conditions of the work- 
piece, leading to its premature failure. 

Much of the work done in the past were to estimate 
the residual stresses experimentally using bending-deflec- 
tion technique or Z-ray diffraction, ¥ith the development 
of numerical techniques like Finite Element Method, attempts 
have been made to predict these stresses theoretically. 

Thermal stresses in an elastic-plastic continuum can 
be obtained from the method suggested by Zienkiwicz et al. 
[18], This method has been applied to evaluate residual 
stresses in grinding by Mishra et al. [17]* In the present 
work, the temperature dependence of Youngs Modulus and 
Poisson ratio [22] are considered and the work has been 
extended to study the effects of depth of cut, table speed, 
and , cooling fluid caa the residual stresses. 
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3.2 MALYSIS 

Most applications of Pinite Element Method to 
Solid Mechanics prohlems rely on the use of variational 
principle to derive the necessary element properties and 
final equations. A displacement field within each element 
is assumed in order to invoke the potential energy principle, 
This is called displacement method or compatihility method 
in Einite Element Aaialysis. 

Consider an elastic body of a given shape, deformed 
by the action of body forces and surface tractions as 
shorn in Eig, 3#1. 

The potential energy of such a body is defined as 
the energy of deformation of the body minus in work done 
on the body by the asternal forces, and is given by 

nCu.v) = 4 // t[6][Bf [C][B3{6 }- 2[6][B]^[C]{ £*}]taA 

A 

- // [E"]{6}tdA- f 6 >dS , 

A 

(3.1) 

where 

t = thickness of the body 

|e®|= coloumn matrix of initial strains, which is 
due to non-uniform temperature distribution 
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, u(x,y) 
^ = i / , 


coloumn matrix of tlie components of displace- 
ment field measured from some datum. 


] = matrix relating strains and displacements (Ref: Appendix) 
j = material stiffness matrix (Ref: Appendix) 

£-1 

J= body force components 

‘]= boundary traction components acting on position 0 of the 
boundary (Pig. 3«l). 


She discrete values of the displacement components 
5 assumed within each element having r nodes and are related 
their nodal values by r interpolating functions Nj^(x,y) as 


o) 


^u(x,y) 





Ik 


I 

i=l 

r 

I 

i=l 



^ = * 

1 

1 — 1 

t — 1 

m 

► 

v^ 

1 i 

1 — « 

1 — 1 

/ 

{ v} , 


( 3 , 2 ) 


)stituting equation (3.2) in equation (3,1) » the potential 
3rgy for an element is 


^(e) 

- // [P*]^®^{6}^®^t^®^dA^®^- {6}^®^ dA^®\ 

^(e) 


( 3 . 3 '! 
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The potential energy of the discrete system assumes 
its minimum value mien the first variation of the functional 
vanishes. At any node q, we have 


3 u 


a 

3 V 


h = {0} = 


// [B]^ [G 1 ^ [B] ^ { 6 t^ dA^ 

(e) ^ Jq ^ q 


A 


// [B]^[C]^®\ e* dA^®^ 

^(e) ^ •‘q ^ o 


ff IS {I'*}^®^t^®^dA^®^ 

ife) 

Cqq 


(3.4) 


where 



coloumn vector of the displacement components 
at node q. 


Equation (3*4) is the Eorce-Displacement relation for the node 
q. In matrix notation, it can he expressed as 
2x2 2x1 2x1 2x1 2x1 2x1 

[K]q{5}q = {J'cJq-*- 


f 


(3.5) 
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where 

[K] 




{i- } 

0 q 


= stiffness matrix at node q 

= If = Initial 

/ e) 

force vector at node q 




2x1 


/jC 


ll 


^■(e) 

vector 


(x,y) { B ®^t^ ^^djL^. = nodal Body force 




C 


// il (xjy) { T® dS ^ = nodal force 

(e) ^ q q 


vector due to surface loading 


{B}g^ = external load vector at node q. 


Equation (3*5) is calculated for all the nodes of an element, 

» 

and than assembled for all the elements to form global stiff- 
ness matrix and load vector* The element chosen here is the 
same triangular element used in the calculation of the tempera- 
tures and the element stiffness matrix is given in Appendix. 
The final equation talces the form 


[Kp] { 6 } = (Fp} 


(3.6) 


4 
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where 


[^3 

{ 6 } 

{^L> 


total stiffness matrix 
displacement vector 
load vector. 


The strains due to temperature rise in the workpiece can be 
written, in case of plane strain as 


r *1 
{ £ } 
^ o 


yo 


■^xyo 


( 3 . 7 ) 


where e and e ^ are the normal strains and Y-v-sm 
xo y 0 xy u 

shear strain. For an isotropic material, the strain vector 
in an element subjected to a temperature rise T is 

[cc All 

{ r= (i+v ) |a ^T r , (3.8) 

where 

V = poisson ratio 

a = coefficient of linear expansion 
AT = temperature increment. 
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The equivalent loads corresponding to the strain vector are 
given by 

J [on "o > Oa'®’ - (3-9) 

The thermal stress at any point in an element, ior elastic 
behaviour of material is obtained from the x-elation 

{S} = [0] [B]{6}[®[ - [0] {c“) , (3.10) 


where 


{S} 









(e) 


element displacements. 


In the present work, Von-mises yield criterion is assumed, 
and the material is assumed to be elastic perfectly plastic. 
For a plane strain case, yield occurs when 

T(T) < i , (3.11) 

where 

Y(T) = yield stress at temperature T 
S = effective stress = /-tC S^+S^-S S +3S^ } 
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longitudinal and transTerse stresses, 
shear stress. 


The check in equation (3«11) is carried out for all the 
elements. In the elements where yield has taken place, the 
excess stresses beyond is calculated and redistributed in 
the following manner : 

(1) Calculate r = « 

S 

(2) Calculate excess stress as (l‘'r)8„^ (l-r)S„_ . 

A y xy 

whore S,, , S, . are the component o of stress for 
x^ y’ xy 

elastic behaviour of material. 

(5) Calculate the corresponding loads 

‘V = - tf, . Pf , (3.12) 

j^{0) 

whore 

{Sq} = vector containing excessive stresses. 

(4) Solve the elastic problem onco again, 

(5) Modify the stresses at all points in the continuum by 
adding the stresses obbainod in step (4) to the previous 
values, 5his process is continued till the diffeir^QO 
botwoon two successive values of stresses reaoh a 

(6) Proceod for all the incremental thermal loads 
from ths ffiOvq?®Qnt of tho wheel, in a similar 
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A computer flow chart depiciting the various steps involved 
in the calculation of residual stresses is given at the end 
of the chapter. 

From the previous chapter, it was found that the 
temperature gradients exist only in a shallow region "below 
the workpiece surface, The displacements at a depth "below 
this affected zone are nogligi"blo compared to those at the 
surface. The displacements on the piano FF (Fig. 2.3) are 
therefore, assumed to "bo zero. Some discretization model is 
usod for "both stress and temperature models [Fig. 2,3]. 

3.3 RSSUITS 

¥ith the temperatures obtained using the theoretical 
analysis in the previous chapter, the stresses in the work- 
piece wore calculated. The input parameters are given in 
Table 2, IB. The results are plotted in Fig, 3.2 - Fig. 3.6. 
The nature of the curve, and tho order of magnitude of the 
stresses match with tho experimental data available [12], 

With tho increase in table speed, residual stresses increase 
(Fig. 3.4), And also with tho increase in depth of cut, 
vhoro is considerable increase in tho peak value of tho resi- 
dual stress (Fig. 3.3). The effect of the cooling fluid on 
the peak values of tensile residual stress is plotted in Fig, 
3,6. The stress history in the worlspiece is shown in Fig.3.2, 
Tho variation of transverse residual stresses axe shown in 
Pig. 5.5. 
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Body 

forc®s 





fig. 3.1 Arbitary two dimensional elastic body 
experiencing surface tractions ond 
body forces. 






0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

Depth below the eurfoce , mm 


Fig. 3.2 Residual stress - history inside the workpiece (set no. 5] 
Table 2.)B. 

1 * 




Strtss , Stress ^ K§/mm* 

Compressive Tensile Compressive Tensile 
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Depth below the surface . jam 


ig. 3.3 Effect of depth of cut on residual stress 
distribution (set no. 1,4.7), table 2.1B 



Depth below the surface , /am 

Fig. 3.4 Effect of table speed on residual stress 
distribution (set no. 1,2,3), table 2.1 



Peak tensile stress , 
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Depth below surface (0.254 mm) 


Fig. 3.7 A typical residual stress distribution. 

After Letner (11^) 

Wheel speed : 30480 mm/sec 
Work speed : 304.8 mm/sec 
Cross feed : 1.27 mm/pass 
Depth of cut '. 0.0508 mm 
Grinding fluid - Oil -water emulsion 
Work piece material ^ 52100 steel 
Hardness : 59 Rockwell ''C” 





i 



0@pth, 0.0254 mm 







0 200 'too 600 JiX) 


Fig.3.9 Assumed mechanical properties of mild steel 
dependent on temperature.[R€ji^ 
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GHAPTER IV 


DISCUSSION OP RESULTS 

Pinite Slomont Analysis yields realistic values, 
provided the material properties and thoir variation with 
temperature are exactly kaowno Pinite Element Analysis 
is based on the concept of discretization, like all other 
numerical approximations. Despite the fact that a solution 
is obtained at a finite number of discrete nodal points, 
the formulation of field varialle model inherently provides 
a solution at all other locations in tho body. The main 
o,dvantago of this method is that complex geometries with 
difficult boundary conditions can bo analysed. But the use 
of this method is restricted to the a,volla,bility of large 
speed computers. 

The present model for obtaining tomperaturo and 
stresses can give realistic values provided tho input heat 
flux Qj,, and its variation along tho contact length is 
exactly known. Prom literature [27], it was observed that 
tho assumption that varies uniformly yields good results. 

It can be soon from tho input data. Table (2, IB), 
that tho forces increase with increase in depth of cut and 
table speed. Tho factors that govern tho forces are : 
individual machine settings and their combinations in various 
forms, dressing techniques, coolants, whoel-worlspieco material 
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combination, rigidity of tlie raachine aiiO so on_, Prom 
equation (2,2) , uitb. ircreatse in force, tns lioat flux Q^, 
increases tnereby increasi'cg the temperatures. This is 
observed in Pigs. 2.4 to 2.14, where the temperatures for 
different depth of cuts -and t.ahle speeds arc plotted, 
Iiicr£'a,siiiQ donth of cut, inci*eases the temperatui'e but with 


the increase in table speed, tlie increase in temperature 
is marginal, lliis is because, wn/bh ijicrease in table speed, 
there is a decrease in the iractir-n of gicindinf;, energy enter- 
ing into the workpiece, which reduces the effect of increase 
in tjie tsuigential force in equation (2,2), Since a is assumed 
to be a constant, the effect of 'oable speed Oii temperature 
is observed to bo ma.rginal. At higher table speeds, however, 
the depth of the affected zone reduces. Pigs. 2,6 to 2,0. 


Since the stresses are resulting from these tempe- 
ratures, they also follow a similar variation with depth of 
cut and table speed. Pigs. 5,3 to 5.4 show the plot of the 


I'esiclual stresses as a function of depth below machined 
surface, for different table sneeds and clenth of cut, They 


aro tensile in nature at the surface and at lower depths they 


tend to become compressive in some cases. At higher depths of 
cut, the peak temperature values obtained are more (Pigs, 2, 8 to 
2.10) due to the fact, P increases. Also a (equation 2.2) 

j p 

increases with depth of cut. So with the incroase in the 
depth of cut, peak values of stress also increase. Prom 
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combiiic/Gioii, of flie raachiiie aiaO so op.. Prom 

« 

oguavion (2,2) , ircroabs in force, ilie licaf flux Qj, 

in crea;' 3 S flLeroby increasing' flio i"oinpera.f iires. flais is 
observed in Pi^s. 2,4 to 2.14, where the temperatures fox' 
different depth of cuts .and table speeds arc plotted, 
Incroasinp, depth of cut, increases the temperature but with 
■ulie increase in table speed, the increase in temperature 
is marginal, lliis is b.^'cause, xrith increase in table speed, 
there is a decrease in. the iractirGi of gzoindi’ng energy enter- 
j.ng into tlie workpiece, which reduces the effect of increase 
in tlie tazigsntial force in enuation (2,2), Since a is assumed 
to be a constant, the effect oi rable speed on temperature 
is observed to bo ma,rgincil« At higher table speeds, however, 
the depth of the affected zone reduces, Figs. 2,6 to 2.8. 


Since the stresses are resulting f3.'om these tempe- 
ratures, they also follow a similar variation with depth of 
cut and table speed, P'igzs, 3*3 to 3*4 sjiow the plot of the 
3 'esiclual stresses as a function of depth below machined 
sincface, ior different table speeds and depth of cut. They 
ai'G tensile in iiauure a.t the surface and at lower depths they 


tend to bscomo compressive in some cases. At higher depths of 

cut, the ]peak temperature values obtained are more (Pigs. 2, 8 to 

2,10) uue to the fact, P increases. A3.so a (equation 2.2) 

P 

increases with depth of cut. So with the inoroase in rhe 
depth of cut, peak va3-ues of stress also increase. 


Prom 
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Any diecrepency In ‘the values can be attribi-'ted to the 
method by which they are obtained* Pinite Element 'jyiethod 
is based on discretization of a body into small regions and 
also, the heat source (wheel) movement is modelled by step 
wise movements. Better results can be obtained by using 
finer Finite Element mesh and smaller movements by the 
wheel. But the computational time will increase. It was 
also assumed that varies uniformly over the contact 
length 1 , but in actual grinding process its variatim is 
unknown. As the cut progresses the grain experiences more 
and more imcut material, and the forces on the grain increa- 
ses, reaching a maximum at the end, when the grain leaves 
the interface surface* So Qj, also varies as the cut progresses. 
Published results [27] indicate that of the three types of 
variations assumed, linear, parabolic, uniform, the last one 
gave fairly accurate results. 



CENTER V 


GOHCLUSIOH ME FUTDRS WORK 

Tlio present model for obtrarxing grinding temperatures 
and residual str esses using Finite Blemont Method yields good 
results. Ihc temperatures tlnat are obtained, are in good 
agreement with the published experimental results. The stress 
distribution that was obtained from these temperatures agrees 
with the results available in literature qualitatively. 

The increase in the depth of cut and table speed 
results in, 

(1) the peak value of the temperature to rise to such an 
extent that thermal damage to the ground surface in 
the form of grinding burn takes place. 

(2) the peak value of the tensile residual stress to 
increase, finally reaching the ultimate strength of 
the material resulting in gi'inding cracks, thereby 
reducing the fatigue strength of the worljpiece. 

At lower table speed the depth of affected zone is increased, 
whereas at lower depth of cut, it decreases. Variation in the 
heat transfer coefficient, so as to simulate the effect of 
coolant, has marginal effect on the peak value of the tensile 
residual stresses. Therefore, for better grinding result^ , 
a lower depth of cut and table speed is preferable. Accuracy 
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of the results can be further improved by employing a finer 
mesli at the ground surface, and decreasing the step length 
of the wheel movement. The present model can be further 
strengthened by including the phase changes that take place 
on the ground sui*face, the heteroginity of the material, 
Further, to obtain realistic values by Finite Element Method, 
the input parameters such as the fraction of the grinding 
energy entering the worljpiece, the variation of heat flux 
along the wheel— workpiece contact zone should be known exactly. 
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APPEISTDIX 


The workpiece is discretized into a finite nninber 
of Constant Strain Triansnlar elements as shown in Pig, 
Consider the figure shown helow, 

y ^ 





where 1,2,3 are the node numbers. The variation of the 
function (temperature or displacement) within this element 
is given by 
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The natural co-ordinates 1^, L 2 » related to the 

cartesian co-ordinates by 
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where 1 - 1»2,3 are the cartesian co-ordinates 

of node i« The area of the element is given hy 

^1 
^2 

The derivatives of the shape functions, which are used in 
assemhling the stiffness matrices are given by 
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The thermal conductivity matrix for an element e is given by 




{k^b| + iCyC|) (k^b2h^+kyC2G^) 

SYlvi ^^x^3 


Thermal capacitance matrix for an el®nent 



is given by 
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Matriic relating the stresses and strains is given by 




he stiffness matrix [Ky,] for an element is given by 
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